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Abstract. Under sufficiently high electric field gradients, 
electron behaviour within exactly perpendicular shocks is un- 
stable to the so-called trajectory instability. We extend previ- 
ous work paying special attention to shortiscale, high ampli- 
tude structures as observed within the electric field profile. 
Via test particle simulations, we show that such structures 
can cause the electron distribution to heat in a manner that 
violates conservation of the first adiabatic invariant. This is 
the case even if the overall shock width is larger than the up- 
stream electron gyroradius. The spatial distance over which 
these structures occur therefore constitutes a new scale length 
relevant to the shock heating problem. Furthermore, we find 
that the spatial location of the short-scale structure is impor- 
tant in determining the total effect of non-adiabatic behaviour 
- a result that has not been previously noted. 

Keywords. Space plasma physics (Shock waves; Numerical 
simulation studies) 



width is larger than the upstream electron gyroradius. This 
behaviour allows an electron to change the kinetic energy as- 
sociated with its gyrovelocity perpendicular to t he magnetic 
field smoothly as it crosses the shock. However. ICold (1 1 9761) 
showed that in the presence of an electric field with constant 
gradient. 
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particles will gyrate at an effective frequency, given by 
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where ^eff and SI are the effective and normal gyrofrequen- 
cies, q is the charge on the particle and m is the particle mass. 
The effective gyrofrequency must then be used in calculating 
the gyroradii, i.e. 



1 Introduction 

Collisionless Shockwaves occur throughout the universe. 
While often cited as the production source for high energy 
cosmic rays, the heating mechanisms that act on the dif- 
ferent sub-populations of particles are still not entirely un- 
derstood. Numerous studies have been conducted into the 
electron heating problem, with the characteristic scale length 
of the shock emerging as an important par ameter govern- 
ing the evolution of the electron distribution fBalik hinet al. . 
[1998; Lembege et al., 1999; Schwartz et al., 201 i|). Addi- 
tionally, despite the amount of work conducted on shock 
scale lengths, there is still a lack of consensus regarding the 
relative scales over which the magnetic and electric fields in 
shocks actually vary. 

Electrons are expected to behave adiabatically, conserving 
their magnetic moments jjLm = W±/B, as long as the shock 
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where r^^^ is a new effective gyroradius and v is the gyrov- 
elocity of the particle. The condition for adiabatic behaviour 
must be revised such that the shock width is much bigger than 
the effective gyroradius. Equation dU shows that, for certain 
values of dEJdx, the effective gyrofrequency can approach 
zero corresponding to an extremely large effective gyroradii. 
Non-adiabatic electron behaviour is therefore possible, even 
at shocks with scale lengths much larger than an upstream 
gyroradius. 

The link between scale l engths and n on-adiabatic heat- 
ing was explored by Balikhi n et al.l (Il993b . The authors con- 
ducted a theoretical analysis of electron trajectories at exactly 
perpendicular shocks and identified the so-called trajectory 
instability. This instability causes two neighbouring electron 
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trajectories to diverge exponentially from each other in phase 
space, causing a breaking of magnetic moment conservation, 
wherever i^L^ < 0, i.e. as long as the following instability 
criterion is obeyed: 



e dE^ 
m dx 



f72>0, 



(4) 



where e is magnitude of the electronic charge, m is the elec- 
tron mass, il is the electron gyrofrequency and dEJdx is the 
electric field gradient along the shock normal. The criterion 
requires that the electric field gradient be above some critical 
value or, equivalently for a given cross-shock potential, that 
the scale length of the electric field be below some critical 
value. This can be fulfilled even if the upstream gyroradius 
is smaller than the shock scale. The authors then showed via 
a series of test-particle simulations that the onset of the tra- 
jectory instabilit y coincides w ith the onset of non-adiabatic 
heating. While Balikhin et alj (11993) draw a strong connec- 
tion between the scale length of the shock and subsequent 
heating, they do not alter the scales of the electric and mag- 
netic fields independently of each other, nor do they study the 
effect of displacing one with respect to the other 

This wor k was subsequently extended into the oblique 
regime by Balikhin et al.l (Il998b . In this paper the authors 
also included terms that account for the changing magnetic 
field, which were previously neglected, and found that the 
divergence in phase space always occurs and that the rate of 
divergence is dependent on the gradients of both the mag- 
netic and electric fields. 



Further relevant work is done by iLembege et al.l (120031) . 
Two approaches were used to analyse the demagnetisation of 
the electrons at the shock front. In the first instance, non- 
stationary and nonuniformity effects were included in the 
form of a fuU-particle self-consistent simulation whilst in the 
second instance these effects have been removed. The au- 
thors found that the fraction of electrons which become de- 
magnetised depends on the nonstationary behaviour found at 
shocks. However, it is difficult to attribute this result to any 
particular process or feature of the shock since it is impossi- 
ble to systematically vary particular variables of interest in a 
full particle code. 

It is clear that the relative scales over which the magnetic 
and electric fields vary have a large impact on the type of 
electron heating that occurs. Indeed, the relative field scales 
of sho cks is a topic w h ich we study within this paper In their 
paper, iBahkhin et al.l (119981) outlined possible relationships 
between the fields, though it is a matter of contention which 
of them occurs in reality, since various observations and sim- 
ulations support differing views. 

It is common that both scales have the same ord er of 
magn it ude in simulations and o bserva ti ons (Balikh in et al. 
19931; iFormisano and Torbertl I l982t iFormisanol Il982 



Balikhin et al.','2002';'Lerov et al 



'l982';'Lieweretal." 



Scholer et al... 2003; Lembege and Dawsoti 




likhin and Gedalin (1994) suggest that the variation of elec- 
tron heating with upstream el ectron thermal Mach number 
Vfiowlvthermai-e, reported bv lSchwartz et al.l (119881) . can be 
recovered in this sim ple configuration. On the other hand, 
Scudder et al.l (Il986l) analysed a shock where the electro- 
static potential varied over a scale larger than the magnetic 
field ramp. 



However. lEselevich et al.l(ll971l) reported on so-called iso- 
magnetic jumps w hich were observed in laboratory plasma 
experiments whilst JHeppner et al.l (Il978l) reported the obser- 
vations from ISEE-1 of large changes in the electric field 
over scales much shorter than the mag netic field ramp . More 
recently. Walker et al. (2004) and Bale and Mozeil (|2007|) 
have shown the existence of short-scale, high-amphtude elec- 
tric field structures or 'spikes' w ithin the overall electric 
field profile with Bale and Mozen (2007) speculating that the 
spikes in the electric field profile may lead to incoherent heat- 
ing of the electrons. 

In this paper, we will show for the first time that this is 
indeed possible. Using test-particle simulations, we will find 
the effect of varying the electric field scale length indepen- 
dently of the magnetic scale length; which has not been done 
before. Additionally, we will vary the location of the electric 
field within the shock. We also investigate the consequences 
of an electric field spike within the shock. In doing so, we 
will demonstrate that these electric field spikes constitute a 
new scale length which is important to the shock heating 
problem, and that its location within the shock layer can dra- 
matically change the amount of heating observed. 

The rest of this paper will be structured as follows. Section 
2 will cover the details behind the simulation, with the results 
and analysis following in section 3. Conclusions follow in 
section 4. 



2 The Simulation 

2.1 Field Profiles 

A test-particle approach, where static electromagnetic fields 
are prescribed, is chosen for t his investigation. The normal- 
isation details can be found in JBalikhin et alJ (1 19931) and are 
briefly reproduced here. Time is normalised to the inverse 
gyrofrequency, Vt^^\ coordinates are normalised to the elec- 
tron inertial length, cuj~^ ; velocity is normalised to the up- 
stream Alfven speed, va', and magnetic fields are normalised 
in terms of the upstream magnetic field strength, B^- The 
field profiles u s ed are based upon the profiles described by 
Balikhin et alj (Il993h . They are idealised versions of ex- 
actly perpendicular coUisionless shocks. The field profiles 
are shown in Fig.[T]and are given by Eq. (JSJ), (|6) and (jT). The 
shock is at rest in the simulation frame, with the upstream- 
pointing normal in the —x direction. 
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Fig. 1. Profiles for the electric, Ex (red curve), and magnetic, Bz 
(blue curve), fields in dimensionless units (see Eq. (|5} - 0). The 
shock width is set by the parameter Db , with the Db = De = 1 
case illustrated in the figure. Note that the E^ scale is negative. 
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Here, A0o is the cross-shock potential and chosen to be 
300eV unless stated otherwise. Ey is constant everywhere 
and calculated from the upstream bulk electron velocity 
and magnetic field strength, Ey — VuBu- We use values of 
Vu = 400kms^^ and i3„ = 5nT which are typical for earth's 
bow shock. When normalised, Ey is equal to the Alfvenic 
Mach number. Ma, which we choose to be Mach 8. De 
and Db are the half-electric and half-magnetic field widths 
normalised to the electron inertial length. Equation (|5) only 
applies within the region of space —De > x > De- Every- 
where outside this region, Er^ = 0. Similarly, Eq. (jT) only 
applies within —Dg > x > Dg, taking the values B^ — I 
for X < —Db and B^—i for x > Db- Adiabatic electron 
behaviour, conserving magnetic moment, would therefore 
correspond to a three-fold increase in the temperature of 
the electron distribution based on the jump in the magnetic 
field. We have chosen to use two scale lengths, De and 
Db, rather tha r i the s ingle parameter, D = De = Db, that 



Balikhin et al.l (119931) use because it is important for this 



2.2 Electron Distribution 

For each simulation run, a Maxwellian distribution at a tem- 
perature of lOeV consisting of 600 electrons is initialised far 
upstream from the shock. Since the shock is exactly perpen- 
dicular, the electrons only require two degrees of freedom in 
velocity space allowing us to set v^ = 0. 

For the purposes of this investigation, the temperature cor- 
responding to the two perpendicular (x, y) degrees of free- 
dom will be defined as follows: 



T = 



2kf 



(v-(v)f). 



(8) 



i.e. the temperature is proportional to the variance of the ve- 
locity vectors of all the electrons in the distribution. In prac- 
tice, the parameter that will be of interest is the heating ratio, 
Rh', that is the ratio of the far downstream electron distribu- 
tion temperature to the far upstream temperature. 

3 Results and Analysis 

To investigate short-scale electric field structures, it will be 
instructive to investigate, separately, the scale and location 
of the cross-shock electric field, E^- We will then move onto 
a final set of simulations in which the cross-shock electric 
field will vary over the same scale as the magnetic field with 
a spike embedded within it to better represent a real shock. 

3.1 Electric Field Scale Length 

For this experiment, we will vary De whilst holding Db 
and the total cross-shock potential, eA0o, fixed. The starting 
shock parameters that will be considered are De — Db — 5 
and eA(/)o = 300cV. This scale length corresponds to a shock 
width, 2Db, of 11.2 upstream gyroradii for a lOeV elec- 
tron. These condition s are adiabatic as shown in Fig. 4 of 
Balikhin et alj (1199 3^ and will be the control case against 



which other simulations are compared. 

Figure |2] shows that as De is decreased, the heat- 
ing remains roughly adiabatic for larger De before in- 
creasing rapidly for scale lengths below De ~ 3. At these 
smaller electric scale lengths, the heating is significantly non- 
adiabatic. By holding the cross-shock potential constant and 
decreasing De, the electric field gradient becomes larger. 
This result should therefore not present much surprise since 
it is already known that the separation of the adiabatic and 
non-adiabatic regimes in perpendicular shocks depends on 
the electric field gradient as given by Eq. (|4]l. According to 
this criterion, the threshold of the trajectory instability occurs 
at De '^3.1 for the parameters of our simulation. 



study that we are able to vary the two scale le ngths indepen- 3.2 Displacement of Electric Field 



dently . These particular forms were chosen bv lBalikhin et al. 



(119931) because they are smooth and well behaved at the 
shock edges and throughout the shock layer. 



Having varied the width of the electric field profile, its po- 
sition relative to the rest of the shock can be altered since 
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Fig. 2. Ratio of downstream to upstream electron temperature as a function of electric field scale length, De- The magnetic field scale length 
is kept fixed at Db = 5. A sketch of the field profiles and their relative scale lengths is shown in the inset. For large De, the heating stays 
adiabatic as the electric field scale decreases. Once the adiabaticity is broken at scale lengths shorter than roughly De ~ 3, however, there is 
a negative correlation between the heating ratio and the electric scale length. 



De is smaller than Db- The inset of Fig. [3] shows the dis- 
placement of the electric field with respect to the magnetic 
field such that their centres of variation no longer coincide. 
For this set of simulations, we fix De = 0.5. As before, 
Db ~ 5 and eA(/)o = 300eV. Figure [3] shows a clear trend 
of higher (lower) heating for displacements towards the up- 
stream (downstream) side of the shock. 

To understand why displacing the electric field would 
change the amount of heating, despite maintaining a constant 
electric field gradient, it is necessary to look at the drifts in 
the system. For the field geometries used, the electrons ex- 
perience an Eyf X Bzi drift in the x direction, together with 
an ExX X Bzi drift and a V|B| drift which are in the +y and 
— y directions, respectively. The Eyf x B^z drift causes the 
electrons to drift through the shock and gain all the poten- 
tial energy associated with the E^ field, i.e. the cross-shock 
potential. This is fixed by the A(/)o parameter. The remain- 
ing two drifts cause the electrons to travel along the shock 
in opposite directions. The V|B| drift is directed such that 
the electrons gain kinetic energy from the motional electric 
field, Ey. Conversely, the E^x x B^i drift is directed such 
that the electrons lose kinetic energy to this field. The lat- 



ter two drifts, in addition to the fixed cross-shock potential, 
therefore determine the net kinetic energy gain of the elec- 
trons as they drift through the shock (IGoodrich and Scudder 
1984). 

It will be useful to compare two limiting cases in our ex- 
planation. The electrons will drift through most of the shock 
before encountering the electric field when it is displaced 
downstream. However, when the field is displaced upstream, 
the electrons will encounter it immediately and gain the en- 
tire cross-shock potential straight away. Since the V|B| drift 
speed is proportional to the kinetic energy of the electron, the 
magnitude of the V|B| drift will be larger in the second case 
as it has gained the energy from crossing the E^ field ear- 
lier. Figure |4] shows the trajectories of three electrons which 
demonstrate this effect. The outer vertical lines represent the 
outer edges of the shock i.e. x — ±Db with the inner verti- 
cal lines representing the edges of the displaced electric field, 
i.e. X = Se^De where Se is the displacement of the electric 
field. All parameters are kept the same with the exception of 
the displacement of the electric field. The electron in panel 
(a) immediately picks up the cross-shock potential energy, 
eAipQ. Initially the E^x x Bzi drift dominates, resulting in 
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Fig. 3. Ratio of downstream to upstream electron temperature as a function of electric field displacement, Se- The electric and magnetic 
field scale lengths are kept at De = 0.5 and Db = 5. A sketch of the field profiles and their relative scale lengths is shown in the inset. 
The displacement of the electric field spike given in terms of Db, i.e. Se = —1 corresponds to the center of variations in the electric field 
coinciding with the upstream edge of the shock layer. The heating ratio is greater for displacements towards the upstream edge of the shock. 
Conversely, when the electric field is displaced towards the downstream end, the heating ratio is lower. 



the loss of some of this energy. The V|B| drift then operates 
in the E^ — region where, due to the enhanced perpendic- 
ular velocity, a large drift velocity results in a net — y drift. 
This corresponds to a large non-adiabatic energy increase. 
Panel (b) is similar but the V|B| drift is less effective since 
the electron spends less time in the post-Ex region, allow- 
ing less time for the V|B| drift to act. In panel (c), there is no 
space for the V|B| drift to act after the electrons have crossed 
the cross-shock potential. The Exit, x B^i drift reduces the 
energy gained from eA0o the most compared to the other 
panels. 

Outside the region in which E^ is non-zero, the energy 
gains associated with the V|B| drift are roughly consistent, as 
expected, with adiabatic compression in the increasing mag- 
netic field. Since this multiplies the existing particle energy, 
it gives the most energy to trajectories suffering early non- 
adiabatic processes as in panel (a). 

To summarise, when the electric field is displaced down- 
stream, the electrons drift through most of the shock adiabat- 
ically, losing energy as a result of the E^x xB^i drift, before 
encountering the non-adiabatic divergence in phase space as 



discussed by iBalikhin et alJ (Il993h . When the electric field 
is displaced upstream, the electrons immediately experience 
the phase space divergence. As the electrons drift through 
the rest of the shock, they will undergo a further expansion 
in phase space due to the V|B| drift, the magnitude of which 
is larger when electric field is displaced upstream. This leads 
to a higher heating of the electron distribution. 

We note that the trajectory instability is an essential in- 
gredient in this non-adiabatic behaviour While the displace- 
ment of the electric field toward the upstream enhances the 
instability by keeping the gyrofrquency, fl, lower in Eq. (4), 
experiments with different values of constant fl (not shown) 
are inconclusive. Thus, we prefer to discuss the non-adiabatic 
behaviour in terms of the various particle drifts. Other ex- 
periments (not shown) in which a field with a larger De is 
displaced remain adiabatic. 

3.3 Shock Spikes 

Whilst it has been instructive to consider these simulations, 
observations show structures with a scale much smaller 
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Fig. 4. Three electron trajectories in the xy plane for different displacements of the electric field profile. In all three cases the magnetic field 
variations occur between the two outer vertical lines. The electric field variations are bound by the two left-most lines in panel (a), the centre 
two lines in panel (b), and the two right-most lines in panel (c). All other parameters are fixed. The drift directions are shown in panel (b). 
The panels show that when the electric field is displaced upstream, i.e. panel (a), the electron will drift in the negative y direction a lot more 
compared to when the displacement is downstream, i.e. panel (c) 



than the total shock width embedded wit hin a larger overall 
electric field profile (IWalker et al.l |2004|) . In one particular 
shock crossing. Walker et al. identified three large-amplitude, 
small-scale structures, the largest of which had a peak mag- 
nitude of around 45mV/m, compared to an average motional 
electric field of around 14mV/m. These were the largest 
amongst the field disturbances observed in the shock and oc- 
curred over the middle 50% of the shock transition. The au- 
thors estimate that the width of these structures to be around 
1 — ScWpg^, with the magnetic field ramp occurring over a 
scale ^10 times this. It shall be the aim of the final set of sim- 
ulations to encapsulate these features, if not the actual values 
themselves. Most importantly, for this particular shock cross- 
ing. Walker et al. report that the structures contribute 40% of 
the total cross-shock potential change. 

The inset of Fig. |5] shows the field profile we used to 
model the electric field spikes. The electric field profile 
shown is constructed by adding together two profiles, both 



described by Eq. (|5]l but with different De values. The im- 
portant features of the profile are De = Db and -D^' '^ <C 
-Dp;. In keeping with th e ratio of the scale lengths observed 



-\spike 



by Walker et all (120041) . D'f'^'' is one tenth of De- While 



IW alker et al.l ( 120041) reported that the overall electric field 
scale is slightly larger than that of the magnetic field ramp, 
the two scales have been kept equal here so as to allow for 
the independent in vestigation of the spik es alone. In any case, 
from the work of Balikhin et al.l (11993 ). we would not ex- 
pect that having De > Db would cause the heating to be 
non-adiabatic, as this would make the electric field gradi- 
ent smaller For simplicity, only one electric field "spike" 
has been modelled. Our base case will again be the adia- 
batic shock where De = Db = 5 and the total cross-shock 
potential is 300eV. We choose 30eV as the cross-spike po- 
tential with the rest of the shock accounting for the remain- 
ing 270eV. We vary the position of the electric field spike 
to investigate its influence on the electron behaviour Figure 
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Fig. 5. Ratio of downstream to upstream electron temperature as a function of electric field spike displacement, S'^J" '^ . The electric and 
magnetic field scale lengths are kept fixed at De = Db = 5 and D^' "^ — 0.5. A sketch of the electric field profile is shown in the inset. 
The magnetic field has been omitted for clarity. The displacement of the electric field spike given in terms of Db i.e. a displacement of -1 
would mean that the variations in the electric field are centered exactly at the upstream edge of the shock. The heating ratio is greater for 
displacement towards the upstream edge of the shock. Conversely, when the electric field spike is displaced towards the downstream end, the 
heating ratio is lower. For all displacements, the heating is non-adiabatic. 



|5] shows that for displacements towards the downstream side 
of the shock, the heating is non-adiabatic, but the amount of 
heating above the adiabatic case is small. For upstream dis- 
placements, there is a much higher non-adiabatic component 
to the heating. 

The conclusions of the previous set of simulations can 
readily be applied here. By embedding a spike into the pro- 
file, there is now a region that satisfies the instability crite- 
rion in Eq. (|4]l when previously there was not, thus pushing 
the shock into the non-adiabatic regime. The same trend is 
noticed, with the heating ratio having high values for dis- 
placements towards the upstream end. However, the heating 
ratio is much smaller in comparison to the previous simu- 
lations; this should not be surprising given that the poten- 
tial drop across the spike is much smaller. Just as before, the 
breaking of adiabaticity occurs earlier for displacements up- 
stream, but the phase space expansion effect due to the V|B| 
drift is not as pronounced since the energy gains associated 
with the spike are smaller. 



lew g^, the width of our spike is at the lim it of 



p^ widths reported bv IWalker et all (120041) . We 



Ai2DE 
the 1 — Sew, 

conclude that in general, the presence of short-scale enhance 
ments to the electric field can push an otherwise adiabatic 
shock into the non-adiabatic regime. The width of the elec- 
tric field spikes therefore constitute a new scale length that 
is important in the study of electron heating at collisionless 
shocks. 

We conducted a final simulation with three spikes at dis- 
placements of -Q.5Db, Q-ODb and 0.5Db embedded within 
an underlying electric field of width 2De = 10. Each spike. 



of width 2D 



spike 
E 



1, contributed 30eV to the cross-shock 



potential with the underlying profile contributing 210eV for 
a total cross-shock potential of 300eV. We find the heating 
ratio for this set-up to be Rh = 4.45, which is not a surpris- 
ing outcome based on our previous results. Figure |5] shows 
that the spikes at Q.QDb and Q.SDb have a minimal effect 
above adiabatic electron behaviour. The non-adiabatic be- 
haviour found here is due predominantly to the upstream- 
displaced spike at —0.5Db- 



V. See, R. F. Cameron & S. J. Schwartz: Electron behaviour due to short-scale structures at shock waves 



4 Conclusions 

It has been the aim of this paper to look at the effect that the 
electric and magnetic field scales have on electron heating at 
collisionless shock waves with a focus on short-scale high- 
amplitude structures in the electric field. Our work builds on 
the existing work of Balikhin et al. (1993) and is motivated 
by the obs ervations of short-sc ale el ectric field structures ob - 
served by 'Walker et alj (l2004 and iBale and Mozej (l2007h . 
Bal ikhin et al. ( 1993) showed that shorter scale lengths can 
lead to incoherent electron heating by satisfying an instability 
criterion with the short-scale electric field spike observations, 
providing a possible means of satisfying this criterion in real- 
ity. We have shown that the presence of small-scale structures 
can indeed push the heating of the electron distribution from 
the adiabatic into the non-adiabatic regime. Specifically, the 
main results of this report can be summarised as follows: 

[1] Shorter-scale electric fields lead to non-adiabatic elec- 
tron behaviour 

[2] The position of these electric fields has been shown, 
for the first time, to play an important role in determining 
the level of non-adiabaticity, with higher non-adiabatic be- 
haviour observed for upstream displacements. This is due to 
the earlier energy gain of the electrons allowing for a large 
magnitude of subsequent V|B| drift. Equivalently, the mag- 
netic moment of the electrons is increased more significantly 
for upstream displacements, allowing the electron to gain en- 
ergy adiabatically in the subsequent magnetic field increase. 

[3] This is true even when considering smaller-amplitude 
spikes embedded within a larger-scale electric field profile, 
provided that the magnitude of the electric fie ld gradient is 



large enough. Such spikes ha ve been observed (IWaUcer et al. 



2004 IBale and Mozed 120071) . 



[4] The existence, scale and location within the shock of 
electric field spikes are therefore important new factors to 
consider in the context of electron heating. 

The next step would be to extend our work into the 
oblique regime. The new scale length associated with the 
electric field spikes is relevant to the discussion of shock 



scale l engths and heating at oblique shocks by Bali khin et al. 



(119981) . Other shock features of interest which could influ- 
ence the electron dynamics include foot and overshoot re- 
gions, as well as the time depen dence of the field profi l es and 
higher-frequency fluc tuations (ILembege and Savoinii 12002 : 



Sundkvist et al.Ll2012h . Finally, it would also be interesting to 



study the electrons within the shock layer, where the strong 
trajectory instability and short scales involved might be ex- 
pected to break the gyrotropy of the distributions. 
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